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In this paper, we report results on the anomalous Hall effect. First, we summarize analyti- 
cal calculations based on the Kubo formalism : explicit expressions for both skew-scattering and 
side-jump are derived and weak-localization corrections are discussed. Next, we present numerical 
calculations of the anomalous Hall resistivity based on the Dirac equation. Qualitative agreement 
with experiments is obtained. 

The anomalous Hall resistivity corresponds to the spontaneous value that takes in magnetic materials the Hall 
resistivity in the absence of applied magnetic field and results from a combination of spin-orbit coupling and spin 

■ polarization. Two different mechanisms responsible for this effect are distinguished: the skew-scattering |jj and the 
side-jump both corresponding to an asymmetric spin-dependent scattering by a potential in the presence of the 

■ spin-orbit coupling. The simplest way to calculate such an effect is to start from the Pauli equation including the 
spin-orbit coupling 
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and to use the Kubo formula for the conductivity 
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where f2 is the volume of the sample, (...} c denotes the configurational average, G + and G~ are the retarded and 
advanced Green's functions G ± (e) = (e ± iO — H) . ef is the Fermi level and v% is the i-component of the velocity 
which contains an additional part due to spin-orbit coupling (the so-called anomalous velocity vso) 



§ : v = £ + v so = P +-^(<rx W). (3) 

*^ ■ The anomalous velocity inserted in (|^) gives the side-jump term whereas the spin-orbit coupling contribution to the 
Green's function (i.e., GqHsoGq where Hso is the spin-orbit coupling and Go the Green's function associated to the 
Hamiltonian in the absence of the spin-orbit coupling) inserted in gives the skew-scattering term provided one 
goes beyond the Born approximation [Q . 

We present a simple application of such a method of calculation in the case of a system with a cubic symmetry and 
an effective magnetic field along the z-axis: thus the anomalous Hall conductivity is equal to the off-diagonal element 
&xy We model the compound in the following way: the total volume of the sample f£ is divided into N cells of volume 
ilo- In each cell, the potential takes a constant value V with a probability distribution P(V) which is characterized 
by its moments (V n ) c = J P(V)V n dV. A suitable choice of the energy origin yields (V) c — 0. We assume that 
there are no correlations in the value of the potential in different cells. In order to achieve analytical calculations, 
we restrict the study to the lowest order with the scattering potential and express the Green's function G in term of 
the average Green's function G in the relaxation time approximation. Neglecting weak-localization corrections, the 
diagonal conductivity is then given by the Einstein relation and for the off-diagonal conductivity, we get 
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for the skew-scattering, and 
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for the side-jump, respectively (detailed calculations are presented in Ref. ||). We have introduced the quantities 

A = h/mc, — fiVp^/Amc 2 . A/j-rn, uj? and aVx\ are respectively, the density of states per unit volume f^o, 

the velocity at Fermi energy and the diagonal conductivity, each for up and down spins. In contrast to a xx and 
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&xy j &xy does not depend on disorder. The Feynmam diagrams associated with these mechanisms are depicted in 
Fig. (1). A simple illustration of these results can be given in the case of a binary alloy A x Bi_ x for which we have 
{V 2 ) c = x(l — x)(ea — £b) 2 and {V 3 ) c = x(l — x)(l — 2x)(sa — £b) 3 where £a(b) is the value that takes the potential on 
side A(B). As a consequence, the anomalous Hall resistivity for skew-scattering is equal to pff ~ a^y I°xx 00 (x — 3x 2 ) 
and for side-jump to pf/ ~ a^y ' /5%. x oc x 2 which is in agreement with the empirical relation pn = ap xx + bp xx but in 
disagreement with the common belief that the quadratic term arises only from the side-jump. 

Within this approach, the weak-localization corrections to the anomalous Hall effect can also be calculated. We 
have considered both Cooperons and Diffusons (see Fig. (1)). The results are the following Q): (i) the Cooperons 
diagrams for the side-jump cancel exactly each other whereas the Diffusons diagrams give a negligible contribution (of 
order (Ti/e f t) A where r is the relaxation time), (ii) the Cooperons diagrams for the skew-scattering give a non-zero 
contribution which includes both spin up and spin down channels. As a consequence, the weak-localization corrections 
to the anomalous Hall resistivity Apn /p° H — Aa xy /<7^ y —2Aa xx /a xx exhibit a strikingly different behavior as compared 
to the normal Hall resistivity for which weak-localization corrections vanish by an exact cancellation of the diagonal 
and off-diagonal parts l|. Due to the presence of two spin channels, such a cancellation can never take place in the 
case of the anomalous Hall resistivity. 

The analytical expressions (Q) and (||) have been obtained in the weak-scattering limit and using the free electron 
approximation. In order to get a more realistic description of the anomalous Hall effect, we have performed numerical 
calculations starting from a tight-binding description and using the coherent potential approximation in order to 
treat disorder. The Green's function is then expressed in term of the t-matrix and the average Green's function: 
G = G + GTG. As a consequence, the conductivity (||) can be split, in the case of a configuration- independent 
velocity, into two different parts 

hi = |^Tr [v t G + (e F )v 3 G-(e F )] + |^Tr [v z G + (e F )T j (e F )G~ (e F )] , (6) 

where Tj is the vertex function equal to (T + G + VjG~T~) c . This formulation is very useful since it allows the exact 
determination of the vertex corrections in the ladder approximation. However, because of the presence of the spin- 
orbit coupling, the velocity is no longer configuration-independent (see Eq. (||)) and it is then not possible to take it 
out of the configuration average (...) c like it is done in Eq. (^). A means to avoid this problem is to start, rather than 
from the Pauli equation, from the Dirac equation 

H = c (a ■ p) + (3mc 2 + V - p B fl (<r • B Gfr ) , (7) 

where a and (3 correspond to the standard Dirac matrices. Indeed, in this description, the velocity operator is simply 
equal to col and then is configuration-independent. 

The system we consider is a ferromagnetic binary alloy A x Bi_ x with an effective magnetic field along the z-axis. 
First, we have expressed Eq. (Q) in the tight-binding approximation for a cubic symmetry, next the self-energy has been 
numerical calculated by the means of an iteration procedure and used to calculate the average Green's function and 
the vertex corrections that we insert in Eq. (^|) in order to get the conductivity tensor. Variations of the anomalous 
Hall resistivity pn — o'xy/o'xx as a function of concentration of disorder and spin-orbit coupling are depicted in 
Fig. (2). We obtain a change of sign of the anomalous Hall resisitivity for a particular value of the concentration 
which corresponds to an exact cancellation of the skew-scattering and the side-jump contributions. Such a change of 
sign has been observed experimentally in PdCo and PdNi alloys || . The variation of the anomalous Hall resistivity 
with spin-orbit coupling is consistent with the Onsager relation since the numerical curve can be very well fitted 
when one consider only odd powers with respect to the spin-orbit coupling. A very good agreement with analytical 
expressions which apply in the weak-scattering limit (close to x = and x = 1) and in the weak-relativistic limit 
(Xso — 0) is obtained. 

To summarize, we have performed both analytical and numerical calculations concerning the anomalous Hall effect. 
Explicit expressions of the skew-scattering and side-jump conductivities have been derived which allows the clarifi- 
cation of the influence of the potential on each mechanisms: whereas the skew-scattering conductivity varies on the 
third order of the potential, the side-jump conductivity does not depend on it. Contrary to what happens for the 
normal Hall resistivity, the weak-localization corrections to the anomalous Hall resistivity do not vanish because of 
the presence of two different spin channels. Numerical results obtained in the case of a ferromagnetic binary alloy 
are in qualitative agreement with measurements in the sense that they show a change of sign of the anomalous Hall 
resistivity. 

V.D. is thankful to J. Barnas for numerous discussions. Part of this work is supported by KBN Grant 
No. 5 P03B 091 20 and NATO Linkage Grant No. 977615. 



2 



[1] J. Smit, Physica 24 (1958) 39. 

[2] L. Berger, Phys. Rev. B 2 (1970) 4559. 

[3] A. Crepieux, P. Bruno, Phys. Rev. B (in press); preprint: ^MKi- ma,t/010137i . 

[4] V.K. Dugaev, A. Crepieux, P. Bruno, preprint: |cond-mat/0103182 , 

[5] H. Fukuyama, J. Phys. Soc. Japan 49 (1980) 644. 

[6] A. Hamzic, S. Senoussi, LA. Campbell, A. Fert, Solid State Comm. 26 (1978) 617. 




FIG. 1. Feynman diagrams for the anomalous Hall effect, (a) corresponds to the skew-scattering, (a') to the side-jump, (b), 
(c) and (d) are the Cooperons corrections to the skew-scattering, (b') and (c') to the side-jump, (d') is the Diffusons correction to 
the side-jump. The curve lines represent the average Green's functions, the dashed lines correspond to the potential (including 
the spin-orbit coupling) and the double straight lines correspond to the ladder part. Symmetrical diagrams have also to be 
considered. 
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FIG. 2. Variation of the anomalous Hall resistivity as a function (a) of the concentration for a spin-orbit coupling Aso = 0.5 
and as a function (b) of the spin-orbit coupling for a concentration x — 0.2. The others parameters are sf /t = 0.2, EA/t — 0.1 
and eb/£ = where t = h/2ma^ (ao is the unit cell parameter) and AeA/t = Aes/t = 0.05 where ±Ae J 4( B j/2 is the value 
that takes the exchange coupling — {iB<J z B c ff on site A(B). The full lines with full circle symbols correspond to the numerical 
calculations and the dashed lines correspond to the analytical results given by Eqs. M) and (ttf). 
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